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Localized electron spins in solid-state systems are widely investigated as potential build- 
ing blocks of quantum devices and computers 1 . In fact, their relatively weak environmental 
coupling makes them ideal candidates for the realization of a long-lived quantum memory 2 . 
However, the spins of spatially separated electrons weakly interact also with each other: 
the implementation of the conditional dynamics thus requires either the use of auxiliary 3 
or "bus" degrees of freedom 4 , or additional resources (e.g., electric fields, more complex 
architectures) that allow a selective enhancement of the effective spin-spin interaction 5 ' 6 . 
Alternative approaches are based on multi-spin encodings of the qubit in coupling-free sub- 
spaces, which render unnecessary the tuneability of the physical couplings during gating 7 . 

The implementation of single- (two-) qubit gates also requires the capability of selec- 
tively addressing single (couples of) spins by means a high spectral and/or spatial resolu- 
tion. Though increased by the qubit encoding in spin clusters rather than single spins 8 , 
the lenghtscales characterizing such local magnetic fields are still highly demanding from an 
experimental point of view. The global manipulation schemes 9 " 11 , where quantum compu- 
tation is performed also in the absence of such individual addressing, thus represent a valid 
alternative, maybe the most likely route in the nearest future. 

In the present paper we propose a possible approach to the implementation of the single- 
and two-qubit gates, based on the qubit encoding in antiferromagnetic (AF) molecular 
rings 12 and on its manipulation by means of either local or global magnetic fields. In 
fact, the characteristic intracluster AF ordering of the spins is shown to allow an inter- 
cluster coupling which is effectively vanishing or finite, depending on the system being in 
its computational space or excited to specific auxiliary states, respectively. As an impor- 
tant point, this switching of the qubit-qubit coupling does not require additional resources, 
nor introduces longer timescales as compared to the ones that are used for performing the 
single-qubit gates. Moreover, the above scheme is shown to provide an efficient means for 
the implementation of a global-manipulation approach, already within a simple ABAB . . . 
linear chain of rings and without requiring a multi-subsystem encoding of the qubit, nor 
the use of a large number of auxiliary units. In both these respects, the existence of excited 
states that comes with the implementation of the qubit in spin clusters rather than single 
s = 1/2 spins represents a crucial resource. The prototypical physical systems we have in 
mind are those substituted Cr-based AF rings characterized by the presence of a ground 
state doublet S = 1/2 which: (i) is energetically well separated from the excited levels, thus 
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allowing to treat the molecule as an effective two-level systems; (ii) possesses a highly or- 
dered, staggered arrangement of the spins, which would result in the effecive cancellation of 
suitably engineered intercluster couplings (see Ref. 13 and references therein). The scheme 
for the quantum-gate implementation, which can in principle be applied to analogous spin 
clusters as well, is first explained in general terms by referring to a model Hamiltonian, 
whose physical justification is given in the final part of the paper. 

Scheme for the quantum-gates implementation. The two-qubit model Hamiltonian 
7i AB = 7i A + 7i B + 7Y^f is defined in the Hilbert space of the effective two-level systems A 
and B, namely {|0), |1) , \i > 1)}a <8> {|0), |1), \j > 1)}b, and reads as follows: 

where the parameters J AB are constant, reflecting the permanent nature of the interclus- 
ter coupling. As schematically shown in Fig. 1(a), 7Y^f leaves unaffected the energy levels 
corresponding to both A and B being in the respective computational spaces, while renor- 
malizing the ones associated to the excited states of B. It should be noted that the above 
coupling is "asymmetric", for the excitation of A to its own higher-lying states does not 
induce any renormalization in the energy levels as far as B is in its computational space: 

<^a<l p>i ^-i ^a>ip<i = ®- (2) 

These Hamiltonian and energy spectrum can be exploited for the implementation of the 
two-qubit gates in a straightforward manner. In fact, the application of a 7r-pulse centred 
on the frequency Tiuj = e^j — ef + J A b B will only excite B from |1) to the excited state \bi > 1) 
if A is in its 1 1) state. A sequence of two such pulses, separated by a time delay r, would thus 
result in a conditional excitation and deexcitation, and more precisely in the application of 
a two-qubit quantum gate defined by the truth table |a) (g) \f3) — > e~ taf36 \a) (g) being 
a : (3 = 0, 1 and 6 = r(e B 1 + J AB — ef )/fi. For 6 = tx the above transformation corresponds 
to a controlled-Z (CZ), that can be combined with two single-qubit rotations of the target 
qubit in order to perform a CNOT 14 . On the other hand, the unconditional dynamics (i.e., 
the single-qubit gates) requires the system to be kept in its computational space throughout 
the time evolution, for the transition energies between the |0) and |1) states of A (B) doesn't 
depend on the setting of B (A) (see Ref. 13 for a detailed discussion concerning the case of 
CryNi molecules) . Therefore, the engineering of a time-independent two-qubit Hamiltonian 
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H such as the shown in Eq. 1, combined with a selective addressing of the spin clusters 
by means of local, pulsed magnetic fields, allows in principle to complete the set of universal 
quantum gates. The local manipulation of the qubits makes it superfluous (though not 
detrimental) for A and B to be physically different from each other. 

The above effective Hamiltonian also allows an efficient implementation of quantum gates 
in the absence of a local control on the qubits. The kind of quantum hardware we will refer 
to in the following consists in a linear array . . . A n --\_B n --\_A n B n An + iB n +-\_ . . . composed by 
two sets of identical molecules (Fig. 1(b)). These are assumed to differ from one another 
with respect to the energy eigenvalues, but to share the qualitative features of the low- 
energy spectrum. Correspondingly, the interaction between A n and its left neighbour B n -\ 
is accounted for by an Hamiltonian 7Y£f which is identical to the 7Y^f of Eq. 1, except 
for the swapping of the molecules' labels. The implementation of quantum gates solely by 
means of global fields requires the use of auxiliary qubits in addition to the logical ones. In 
the following, the linear array will be assumed to consist in an alternated sequence of the 
two: the latter encode the quantum information, whereas the former are all set to |0), with 
the exception of the so-called "control unit" (CU), which is set to |1} 9 . As demonstrated in 
Ref. 11, the capabilities required for the gates implementation with such a quantum hardware 
consist in: (i) the separate switching, on and off, of all the A n <-> B n and B n -\ A n 
couplings, and the application of Ucnot to the corresponding set of coupled cells; (ii) the 
performing of general rotations in each of the two subsets A and B. More specifically, the 
single-qubit gates can be performed with an overall number of qubits which is twice that of 
the logical ones, whereas the complexity of the gating and the number of auxiliary qubits 
remarkably increase as far as the CNOT is concerned 15 . 

In the present case, the difference between the transition energies of A and B and the 
directionality of the effective coupling (Eq. 2) allow both the above conditions to be met. 
The conditional excitation of the molecules A (B) to their auxiliary state \a{) (\b\)) can 
be exploited for the implementation of the Ucz involving each A n (B n ) and its left neigh- 
bour B n _i (An), as shown in the previous paragraph (requirement (i)). For the interaction 
Hamiltonians and Tif^ to be effectively turned off and the single-qubit rotations to 
be performed, instead, it suffices to manipulate the spin clusters within the respective com- 
putational spaces (requirement (ii)). Moreover, the existence of a futher excited state in at 
least one of the molecules (e.g., \b 2 > 1) in B) avoids the increase of additional auxiliary 
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qubits for the CNOT implementation. In practice, the sequences of the SWAP gates that 
move the information within the array rigidly shift the states of the identical cells; being 
the logical qubits all encoded alternatively in A or in B and the physical interaction only 
effective between nearest neighbours, the main obstacle in performing the CNOT essentially 
consists in bringing the control- and the target-qubit next to each other. The specific re- 
source that is provided by the present spin clusters is the above mentioned state \b 2 ): our 
aim is that of using \b 2 ) in order to block, and eventually invert, the motion of the qubits it 
collides with, thus provoking two or more logical qubits to be encoded on adjacent cells. In 
fact, the Uswap rotation can be decomposed into a sequence of three Ucnot 14 '- 

TjA n B n TjA n B n TjB n A n TjA n B n fn\ 

U SWAP — u CNOT u CNOT u CNOT-> \°) 

where A n (B n ) acts as a c-qubit (t-qubit) in the first and third Ucnot-, and as a t-qubit 
(c-qubit) in the second one. A Ucnot can in turn consist in the combination of a controlled- 
Z (JJcz) and two rotations of the t-qubit, as shown in Fig. 2(a). As a key point, the 
initial setting of B n to its \b 2 ) state renders uneffective the transformations (i.e., puts out of 
resonance the transitions) enclosed in the dotted boxes; it is easy to verify that the remaining 
ones result in an identical transformation. In the same way, the application of U^r^p™ does 
not change the state of Ai-i (nor that of B n ) if B n is initially set to \b 2 ). 

In Fig. 2(b) we show how to exploit this block for the implementation of a CNOT, where 
\4> no ) and \4>no-i) act as a c- and a t-qubit, respectively. The required steps are the following: 
(Si) after the CU has been positioned on the right-hand side of the c-qubit, a transformation 
U 2 = 1 1X^2 1 + \b 2 ) (1 1 is applied (7r-pulse with central frequency Two = — ef ). As a result, 
the only B cell that is excited to the auxiliary state \b 2 ) is the CU B no , whereas all the 
others remain in the |0) state 16 ; (S 2 ) a rotation Ug{^ P swaps all the qubits but A no and 
£> no , thus allowing the collision of the control and target (logical) qubits, |0 no ) and |0 no _i); 
a transformation (S 3 ) ^ctvot 1 1S performed, which is only effective on |</> no _i), being all 
the cells A n -^ no in their |0) state; (S4) a second Ug W B ^ P , as in S 2 , brings all the qubits 
back to their original positions; (S5) a second U 2 = |1)(6 2 | + |&2)(1|> as in Si, reinitializes 
the CU to its |1) state. 

Physical implementation of the qubit. We now consider in greater detail the physical 
implementation of the qubit by referring to a class of prototypical systems, namely the 
substituted Cr^Ni rings (with odd x), where the presence of the Ni ion typically results in 
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the formation of a ground-state doublet. The precise identification of a spin cluster that 
meets all the above requirements is beyond the scope of the present paragraph, which is 
rather aimed at demonstrating the suitability of these specific molecules at a qualitative 
level, and at pointing out possible solutions for future molecular engineering. The single- 
qubit terms, which we have so far denoted by 7i x=A,B , correspond to the following spin 
Hamiltonian of the related Cr^Ni molecule 12 : 

x+1 x+1 

i=l i=l 

x+1 x+1 

+ ^d i [(s z ,if - s,{s, + l)/3] + ^ s * • D *i • s i> ( 4 ) 

i=l i<j=l 

J i=2 - x +i = J and J XjX+ i = J', g i= i- x = g and g x+1 = g', s i=1 - x = 3/2 and s x+1 = 1, being 
i = x + 1 the site of the Ni substitution. As discussed in detail in Ref. 13 for the case 
x = 7, the ground-state doublet that forms the computational basis {|0}, |1)} approximately 
corresponds to {\S = 1/2, S z = —1/2), \S = 1/2, S z = +1/2)}; these states are characterized 
by an antiparallel alignment of the neighbouring spins, which we propose to exploit in order 
to obtain the effective cancellation of the intercluster coupling. The auxiliary states |&i) 
and 1 62) are respectively identified with the first excited states M = —3/2 and M = —1/2 
(bi,2 = 2, 3), belonging to the multiplet S = 3/2. 

The realization of an intercluster coupling has been demonstrated in similar systems 17 
(ferromagnetic spin clusters), as well as quantum- mechanical entanglement 18 . For the sake 
of clarity, we shall assume in the following that an effective coupling between the rings A 
and B may result from one or more links that connect individual spins belonging to the 
two of them. An Ising-like interaction between spin i of cluster A and spin j of cluster B, 
namely 7Y^f = J^s^sfj, results in an analogous coupling between A and B: 



AB = 4 B (<x\\sU\<*)msi0) S *A 

a!i ' VSa(S a + 1)(2£U + 1)S B (S B + WSb + 1) ' 



where S^ a (Sfp) is the average value of Sf(Sf) corresponding to state \a) of A (\/3) of B), 
and (0||s^ B ||0) = (l||s^ s ||l), being the reduced matrix elements constant within each S- 
multiplet. In presence of more than one spin-spin interaction, Eq. (5) can be generalized by 
summing over the i and j indexes. In particular, we will consider the case where the spin i of 
cluster A is linked to two adjacent spins, j and j + 1, of B, with J^ B = J A f+\ (see Fig. 1(c)). 



TABLE I: Reduced matrix elements (/3\\s Zik \\/3) in the Cr 5 Ni (up), Cr 7 Ni (middle) and CrgNi 
(down) rings. The index values k = 6 (up), k=8 (middle), and k = 10 (down) correspond to the 
Ni substitution. 



p 


S z ,l 


Sz,2 S Zj3 


s z,4 s z,5 


Sz,6 &z,7 


S z ,8 S Z: g 


s z ,io 


0,1 


1.48 


-1.17 1.40 


-1.17 1.48 


-0.80 - 






bi,b 2 


1.51 


-0.13 1.28 


-0.13 1.51 


-0.17 - 






0,1 


1.39 


-1.13 1.30 


-1.12 1.30 


-1.13 1.39 


-0.77 - 




hM 


1.40 


-0.30 1.09 


-0.11 1.09 


-0.30 1.40 


-0.39 - 




0,1 


1.33 


-1.10 1.23 


-1.08 1.21 


-1.08 1.23 


-1.10 1.33 


-0.75 


h,b 2 


1.31 


-0.44 1.01 


-0.21 0.90 


-0.21 1.01 


-0.44 1.31 


-0.38 



As a consequence of the approximately staggered (AF) arrangement of the reduced matrix 
elements characterizing \/3 = 0, 1) (see Table I and relative discussion), the sum of the two 
individual couplings results in a cancellation of the effective one between the rings: 

WWsfWP) = -(P\\sf +1 \\P) ^ J$* = far = 0,1. (6) 

On the other hand, the excitation of B to a state lacking such intracluster AF arrangment 
switches on the effective coupling to A (see Fig. 1(d)). It should be noted that the key 
result, namely Eq. 6, doesn't depend neither on the assumed Ising-like form of W^f , nor 
on the presence of two disinct links: an Heisenberg coupling or a single link symmetrically 
coupled to ■ and sfj +1 would lead to the same result. 

In order to verify to which extent Eq. 6 holds in realistic Cr^Ni rings, we have performed 
detailed calculations for x = 5, 7, 9 (Table I). While a relevant increase is actually achieved 
by exciting the system from |0) or |1) to the vanishing of J^p in the computational basis 
would require slightly asymmetric links, such that Jff I Jtf+\ = — (/^ll s j+i||/^)A/^ll s jll/^) = 
\(P\\s j+1 \\P)/{P\\s j \\P)\(P = : l). 

As already mentioned, a crucial role within our scheme is played by the (excited) auxiliary 
states of the molecules. The same possibility of inducing direct transitions between these 
and the ground-state doublet by means of an oscillating magnetic field is a non-trivial issue, 
for it requires non-vanishing matrix elements (f3\ Ym=i 9i s i\P')j with f3 = 0, 1 and (5' — bi, b 2 - 
In fact, in Cr 7 Ni these are approximately proportional to gm — Qcr — 2.2 — 1.98, and are thus 
finite only due to the presence of the Ni substitution; besides, they depend on the value of the 
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static magnetic field B (for B = 2 T 13 , for example, |(0| ^ i=1 giS Xji \l) / {0\ ]T. =1 9iSx,i\ b i) I = 
3.8). From these calculations, and from others performed on analogous systems (e.g., Cr 7 Fe 
and Ni 7 Cu), we conclude that, as far as the impurity-Cr exchange is comparable to the 
Cr-Cr one, the main physical properties required by the discussed approach are satisfied. 
Therefore, a possible implementation of the linear A-B array 19 may be composed by Cr x Ni 
rings with two different values of Cr ions (xa ^ x B ), or by two among the existing variants 
a given of Cr x Ni (xa = x B ), characterized by different values of the physical parameters J 
and g. 

In conclusion, we have proposed a scheme for the implementation of quantum gates in AF 
spin clusters which does not require the intercluster interaction to be tuneable during gating. 
Having in mind the prototypical case of the Cr-based rings, we have shown how specific 
auxiliary states allow in principle to control the effective coupling between the qubits, and 
to implement global-manipulation approaches, with a minimum number of auxiliary qubits, 
smaller than the one required by single-spin encodings. 
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FIG. 1: (a) Schematics of the lowest energy levels of the AB system with and without (left- and 
right-hand sides, respectively) the effective intercluster coupling. The levels corresponding to the 
computational space do not undergo energy renormalizations; being S^ = —S^i, JqE = —J\b x 
(see Eq. 5). (b) Sketch of the linear array of rings ABAB ... , with asymmetric couplings between 
each two neighbouring units, (c) Intracluster ordering of the electron spins characterizing the 
ground-state doublet of a prototypical AF ring, namely the substituted CryNi. Possible intercluster 
molecular links (solid lines), locally coupling spin i of ring A with j and j '• + 1 of B. (d) The excited 
state \b±) of B lacks the above antiparallel arrangement of the adjacent spins. 
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FIG. 2: (a) Decomposition of Ugy^ P in terms of Ucz and single-cell rotations (Y and Z correspond 
to the Pauli matrixes a y and a z , respectively). The rotations enclosed in the dotted boxes are 
uneffective if B n is initially set outside from its computational space, to the auxiliary state 1 62) - 
The resulting unitary transformation can thus be written as Uqwap = |00)(00| + |11)(11| + |01>(10| + 
|10)<01| + (|0)<0| + |1)<1|)®|&2><&2| • (b) CNOT gate applied in five steps (Si_ 5 ) to \</) no ) and \</) no -i), 
acting as control- and target-qubit, respectively. On top of the wires we report the quantum states 
of the cells after each step (for simplicity, such notation refers to the case where the subsystems 
are never entangled with each other; we can thus write \4> m ) <8> Wno-i) — UcNOT\4>m) ® Wna-l))- 
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